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Dynamical dark matter (DDM) is a new framework for dark-matter physics in which the dark sec- 
tor comprises an ensemble of individual component fields which collectively conspire to act in ways 
that transcend those normally associated with dark matter. Because of its non-trivial structure, 
this DDM ensemble — unlike most traditional dark-matter candidates — cannot be characterized 
in terms of a single mass, decay width, or set of scattering cross-sections, but must instead be 
, described by parameters which describe the collective behavior of its constituents. Likewise, the 

components of such an ensemble need not be stable so long as lifetimes are balanced against cosmo- 
logical abundances across the ensemble as a whole. In this paper, we investigate the prospects for 
£N) \ identifying a DDM ensemble at the LHC and for distinguishing such a dark-matter candidate from 

» | ■ the candidates characteristic of traditional dark-matter models. In particular, we focus on DDM 

O |' scenarios in which the component fields of the ensemble are produced at colliders alongside some 

' number of Standard-Model particles via the decays of additional heavy fields. The invariant-mass 

distributions of these Standard- Model particles turn out to possess several unique features that can- 
not be replicated in most traditional dark-matter models. We demonstrate that in many situations 
it is possible to differentiate between a DDM ensemble and a traditional dark-matter candidate on 
the basis of such distributions. Moreover, many of our results also apply more generally to a vari- 
ety of other extensions of the Standard Model which involve multiple stable or metastable neutral 
particles. 



I. INTRODUCTION 



Recently, a new framework for dark-matter physics has been proposed P, Q . This new framework is called "dy- 
namical dark matter" (DDM), and previous discussions of DDM have focused on its overall theoretical properties [l[ 
as well as on the theoretical and phenomenological implications of a specific model 0, H| that was constructed within 
this new framework. In this paper, we shall return to consideration of the general DDM framework as a whole, and 
discuss the detection prospects for a broad class of DDM models at colliders, and in particular at the Large Hadron 
Collider (LHC). 

As discussed in Refs. [H, HJ, the central hallmark of the dynamical dark- matter framework is that the dark sector 
consists not of one (or merely a few) stable dark-matter particles, but rather an ensemble of constituents which act 
£N) ■ collectively in ways that transcend the physics normally associated with more traditional dark sectors. Such a DDM 
ensemble is not randomly assembled, but instead has certain internal structures which guarantee its phenomenological 
viability. For example, stability on cosmological time scales — normally considered to be a sacrosanct property for 
traditional dark-matter candidates — is not a requirement for such an ensemble. Instead, lifetimes are balanced against 
cosmological abundances across the different constituents of this ensemble in such a way that those components with 
larger decay widths (and consequently shorter lifetimes) necessarily have smaller cosmological abundances, and vice 
versa. Indeed, as discussed in Refs. @, Q, this balancing represents a novel way of satisfying phenomenological bounds 
on the dark sector without imposing stability as a whole, and ultimately represents the most general dark sector that 
can be imagined. 

Most traditional dark-matter candidates can be characterized in terms of their masses and couplings to Standard- 
Model (SM) states, and indeed most phenomenological bounds in the dark-matter literature are phrased in terms 
of constraints on these variables Q. By contrast, in the DDM framework the dark- matter "candidate" is the entire 
ensemble, and the parameters which ultimately characterize a DDM ensemble describe not only the couplings of its 
individual constituents to SM states, but also the internal structure of the ensemble itself. In general, such an internal 
structure might consist of relationships between the masses, relic abundances, and couplings of its components. As a 
result, the natural parameters which characterize DDM models and their phenomenology are fundamentally different 
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from those which suffice to describe traditional dark-matter candidates. Indeed, the most fundamental parameters 
which characterize a DDM ensemble are those which describe how quantities such as the constituent-particle masses, 
abundances, decay widths, and cross-sections scale with respect to one another across the ensemble as a whole. This 
clearly represents a new way of thinking about a dark sector, but the need for such an approach is one of the primary 
features of the DDM framework. 

In this paper, we shall focus on the implications of these fundamental differences for the collider phenomenology 
of DDM ensembles, and in particular for their discovery potential at the LHC. The canonical channels in which one 
generally expects to obtain evidence for a DDM ensemble involve substantial missing transverse energy (hereafter 
denoted fx) — just as is the case for a traditional dark-matter candidate. It is therefore crucial to develop strategies 
for distinguishing between these two classes of models once an excess in one or more of these fx channels has been 
identified. One such strategy, which was discussed in Ref. |3j, is to search for correlations between fx signatures 
and signals in other channels to which the DDM ensemble might simultaneously give rise — channels not normally 
associated with dark matter. Indeed, in DDM models, only those portions of the DDM ensemble which are stable on 
collider time scales contribute to missing-energy signals at colliders. By contrast, it is possible that other portions of 
the DDM ensemble will have much shorter lifetimes and therefore manifest themselves in different channels entirely — 
channels that may not involve fx whatsoever. As discussed in Ref. 3j, a wide variety of DDM scenarios generically 
give rise to observable excesses in both classes of channels simultaneously. However, in a variety of other DDM 
contexts, such correlations among channels may not be possible due to the particulars of the model and the lifetimes 
of the constituent fields in the ensemble. In such contexts — and especially when those constituent fields manifest 
themselves only through fx signatures — alternative strategies for identifying DDM ensembles are necessary. 

Such alternative strategies will be the primary focus of this paper. In particular, one of our central aims is to 
demonstrate that in many cases it is possible to distinguish DDM ensembles from more traditional dark-matter 
candidates based solely on results from channels in which the dark-matter particles manifest themselves as fx alone. 
For example, the invariant-mass distributions of SM states produced in association with the constituent fields in a DDM 
ensemble by the decays of heavy particles can exhibit qualitative features which transcend those usually associated 
with traditional dark-matter candidates. As we shall see, such invariant-mass distributions can therefore provide a 
powerful experimental discriminant between DDM ensembles and these traditional candidates. Furthermore, because 
this technique relies solely on signatures in detection channels in which the component fields of the DDM ensemble 
appear as fx, these signatures are insensitive to the precise lifetimes of any component fields which leave an imprint 
on the invariant-mass distribution in question. Hence, they are likewise insensitive to the characteristic instability of 
the DDM ensemble as a whole, and this remains true provided that these fields are all sufficiently long-lived so as 
not to decay within the detector volume. These aspects of our discussion therefore have a broad applicability even 
beyond the context of the DDM framework, and apply quite generally to any multi-component dark-matter scenario 
involving an additional heavy "parent" particle which decays to final states involving the dark-sector fields, or to any 
scenario involving multiple metastable neutral fields and such a parent particle. 

In order to explicitly illustrate these themes, in this paper we shall consider the case in which each parent particle 
decays directly to a single constituent field within the DDM ensemble along with a pair of strongly-interacting SM 
particles (i.e., quarks or gluons). Even in this simplest non-trivial case, we shall demonstrate that there exist a range 
of characteristic features which are imprinted on the invariant-mass distribution of the two resulting jets and which 
can permit one to distinguish a DDM ensemble from any traditional dark-matter candidate. Moreover, we shall show 
that there exist production mechanisms for the parent particles (such as pair-production via strong interactions) with 
event rates sufficient to enable such a differentiation at the 5er significance level within the first 30 fb _1 of integrated 
luminosity at the sfs = 14 TeV LHC. 

This paper is organized as follows. In Sect. [HI we outline the general aspects of the collider phenomenology of DDM 
ensembles and discuss the strategies for indirectly observing those ensembles in different classes of DDM models from 
the perspective of effective-operator analysis. In doing so, we devote particular attention to models in which the 
component particles in the DDM ensemble are produced via the decays of a heavy parent particle. In Sect. IIII1 
we calculate the invariant-mass distributions associated with pairs of SM states produced in conjunction with each 
dark-matter particle via this mechanism and compare them to the invariant-mass distributions obtained in theories 
involving only a single particle stable on collider time scales. In Sect. IIV1 we assess the statistical significance with 
which such non-traditional invariant-mass distributions can be distinguished from those in more traditional dark- 
matter models. In Sect |Vj we provide an example of one production mechanism which naturally provides event rates 
of the order required for such a differentiation within the first 30 fb _1 of integrated luminosity at the y/s — 14 TeV 
LHC — namely, the pair-production of strongly-interacting parent particles with masses near the TeV scale. Finally, 
in Sect I VII we provide an assessment of how various subtleties associated with certain production mechanisms for 
the parent particle (such as the combinatorial background associated with incorrect pairings of final-state jets for 
processes which yield more than one such parent particle per event) are expected to impact our results. We also 
comment on the broader applicability of our results to scenarios outside the DDM framework. 
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II. DYNAMICAL DARK MATTER AT THE LHC: GENERAL CONSIDERATIONS 

One of the hallmarks of the DDM framework is that the dark sector consists of an ensemble of particles Xn, where 
n = {0, . . . , N}, with N presumed to be relatively large, i.e., N ^> 1. For convenience, we shall label these particles 
in order of increasing mass, i.e., m n+1 > m n . Moreover, as discussed in Refs. [H-Q, these particles Xn exhibit a broad 
spectrum of decay widths T n which scale inversely with their corresponding cosmological abundances tt n . As a result 
of these different decay widths, different Xn may manifest themselves in qualitatively different ways at colliders — 
even in situations in which the Xn all have similar quantum numbers and are therefore produced via similar processes. 
Those Xn with lifetimes r„ > 10~ 10 s are stable on collider time scales and appear in a collider detector as |?t- By 
contrast, any of the Xn with lifetimes r„ < 10 -10 s decay within the detector volume. Provided these rapidly-decaying 
Xn decay predominately to final states involving SM particles, evidence for a DDM ensemble could potentially be 
obtained via the observation of signals in complementary channels which individually provide evidence for Xn within 
either of these r„ regimes. Such signatures are discussed in Ref. 0]. However, for models in which t„ > 10 -10 s for 
all Xn — or for models in which the Xn that decay within the detector volume are not collectively produced at rates 
sufficient to yield observable effects at colliders — such multi-channel correlation techniques are not useful. In such 
cases, alternative strategies must be found for distinguishing DDM models from more traditional dark-matter models. 

In this paper, we will examine one such strategy, which is applicable to a broad class of DDM models possessing 
two key characteristics. First, in addition to the constituent fields Xn of the DDM ensemble, the field content of the 
model must include one or more heavy particles ip which can be produced at a substantial rate at a hadron collider. 
For example, if the tp transform non-trivially under the SM SU(3) C gauge group, they can be produced copiously 
via their interactions with the quark and gluon fields. Second, these additional, heavy particles must decay with a 
sizeable branching fraction into final states including at least two SM fields, along with one or more of the Xn- Decay 
topologies of this sort arise generically, for example, in a specific class of DDM scenarios in which both the parent 
particle ip and the constituents of the DDM ensemble are charged under an approximate symmetry. However, in such 
cases it is an important property of the general DDM framework that such a symmetry is neither required nor need 
be preserved exactly. Therefore, the ensemble constituents need not ultimately be stable. 

We shall demonstrate that it is possible to differentiate between traditional, single-particle models of dark matter 
and multi-component scenarios, such as those which arise in the DDM framework, by examining the invariant-mass 
distributions of the SM fields produced by ip decays. Of course, in cases in which only a small number of the Xn are 
kinematically accessible in if) decays, these invariant-mass distributions are distinguished by the presence of multiple 
kinematic edges. Note that similar features arise in other contexts as well, most notably that in which a parent 
particle can decay into final states involving different multiplicities of the same stable dark-matter particle [U, Q . By 
contrast, in cases in which the number of kincmatically-accessible Xn is large and the decay phenomenology of the %j) 
particles depends more sensitively on the full structure of the DDM ensemble, qualitatively different features emerge. 
In particular, while individual kinematic edges are no longer manifest, the invariant-mass distributions can exhibit 
distinctive shapes not realized in single-particle dark-matter scenarios. These distributions can therefore provide a 
powerful experimental discriminant between DDM ensembles and more traditional dark-matter candidates. 

We note that the technique described above has a broad range of applicability because it is not predicated on 
the observation of signals of both collider-stable and promptly-decaying Xn, but rather of collider-stable fields alone. 
Moreover, for the same reason, this technique can also be applied more broadly to a wide variety of multi-component 
dark-matter models, or to other scenarios which involve multiple particles which are stable on collider time scales. 
Such situations can arise in certain limits of traditional dark-matter scenarios in which the dark-matter candidate is 
stabilized by a parity symmetry such as i?-parity in supersymmetric models or KK parity in higher-dimensional 
theories in which the SM propagates in the bulk ja4l0l]. For example, while only the lightest parity-odd particle 
is absolutely stable in such scenarios, situations can arise in which the decay rates of heavier parity-odd particles 
with similar quantum numbers are suppressed (either by kinematics or by some additional consideration) to such an 
extent that they are also stable on collider time scales. In such cases, the invariant-mass distributions associated 
with the decays of even heavier fields in the theory include contributions from final states involving all such stable or 
metastable particles. These distributions can therefore yield valuable information about the overall coupling structure 
and mass spectrum of the theory. 

In order to analyze the collider phenomenology of DDM ensembles, it is first necessary to characterize the properties 
of such dark-matter candidates in a straightforward and physically meaningful manner. Indeed, a DDM ensemble is 
not a single particle or a small group of particles, but rather a vast collection of individual states whose collective 
properties dictate the dark-matter phenomenology of the model. Taken together, therefore, an ensemble of such 
states constitutes a dark-matter candidate which cannot be characterized in terms of a single well-defined mass, decay 
width, or set of cross-sections for processes involving SM fields. Instead, a DDM ensemble is more aptly characterized 
by parameters which describe its aggregate internal structure. Such parameters may include, for example, scaling 
exponents in certain relations between masses and cross-sections, or between cosmological abundances and decay 
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widths, that hold across the DDM ensemble as a whole. More specifically, one natural set of such parameters might 
include the density of states expressed as a function of the masses m„ of the x„, and a set of exponents which describe 
the scaling behavior of the couplings of the \n to other fields present in the theory. 

In any arbitrary DDM model, such couplings can be described by a set of operators oi™j„ 3l ,„, where the indices rii 
indicate the dark-sector fields Xn, involved and where the index a labels the operator in question. We shall let d a 
denote the mass dimension of the operator oi"] n2 ,...- We shall assume for the sake of simplicity in what follows that 
the only 0i"| nS) .„ which play a meaningful role in the collider phenomenology of the DDM ensemble are members of 
the subset O n a ^ of operators which involve only a single Xn- We therefore have an effective interaction Lagrangian of 
the form 

N 
a n— 

where A is the cutoff scale of the effective theory, and where c na is the dimensionless operator coefficient associated 
with . In addition, we shall also assume that all the constituent particles in the ensemble have the same quantum 
numbers, and thus that the set of O n a ^ consistent with the symmetries of the theory shares a common operator 
structure for all n. In other words, for any given a, the only n-dependence appears in the mass m n of the field Xn and 
the coefficient c na . Furthermore, we assume that the distribution of c na across the ensemble depends solely on m n , 
and that both this distribution and the mass spectrum of the ensemble exhibit general scaling relations of the form 

m n = mo + n s Am 



Cna = c 0q , (2.2) 

V m o/ 

where the (positive) mass-splitting parameter Am and the scaling exponents 5 and j a are free parameters. Note that 
5 > by construction. All of these assumptions can, of course, be relaxed; however, the qualitative results obtained 
below still continue to hold in more general cases as well. Note that the O n °^ are ultimately responsible both for the 
production of the Xn at colliders, and for their decay phenomenology — both in terms of their lifetimes r„ and in 

(a) 

terms of their decay products. Of course the O n which play a dominant role in the production of the Xn need not 
be the same as those which play a dominant role in determining their decay properties. 

In this paper we shall focus primarily on the Xn which are stable on collider time scales, as discussed above; 
hence our primary concern in what follows will be with production rather than decay. In fact, given the scaling 
relations and parameterizations in Eqs. (|2.1[) and (I2.2[) . we can now characterize precisely which components of a 
given DDM ensemble contribute to fx signals once produced. Assuming that a given Xn decays primarily to final 
states comprising SM fields the sum of whose masses is much smaller than m n , we find that r„ scales roughly like 

m 2J a -7 



E 



As discussed above, the distribution of T n across the DDM ensemble has a significant impact not only on the cos- 
mological aspects of that ensemble, but on its collider phenomenology as well. When the lifetime t„ = l/T n of the 
particle is short on collider time scales (r„ < 10~ 12 s), the particle decays promptly within the detector. When the 
lifetime is long (r„ > 10~ 10 s), the particle decays outside the detector and appears as fx- I n the intermediate region 
(10 -12 s < r n < 10~ 10 s), the particle may give rise to a displaced vertex. Thus Xn within different ranges of m n in 
the ensemble may manifest themselves either as fx or via their decay products, depending on the value of r„. For 
example, in the case in which a single operator On governs the decay width of all relevant Xn, Eqs. (|2.2|) and (|2.3j) 
together imply that the requirement for appearing as fx is 

T 1 + n s < 6.58 x 1CT 15 GeV , (2.4) 



m / 

where 7 a here denotes the particular scaling exponent associated with the operator which effectively controls the 
decay width of %„, and Tq is the decay width of the lightest state in the ensemble. 



III. IMPRINTS OF DDM ENSEMBLES IN KINEMATIC DISTRIBUTIONS 



By and large, when a dark-sector field is produced at the LHC via the decay of a heavy parent particle tjj, multiple 
SM fields can also be produced via the same decay. One useful variable that can assist in distinguishing models 
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involving a single dark-matter candidate from those involving more complicated dark sectors in scenarios of this sort 
is the invariant mass of the additional SM particles produced by decays of the parent particle. The identification 
of features in invariant-mass distributions has been shown to be effective in differentiating between the different 
symmetries which might stabilize a dark-matter particle in traditional dark- matter models [5|. In this paper, we show 
that it is also effective in distinguishing between traditional and DDM dark sectors. 

For purposes of illustration, we focus on the case in which each ip decay yields two SM particles and a single Xn- 
Indeed, this is the simplest case in which a non-trivial invariant-mass distribution is obtained for the SM particles 
from the decay of a single ip. Furthermore, since strongly-interacting particles can be produced copiously at colliders, 
we shall focus on the case in which ip carries SU (3) c charge and decays to a pair of SM quarks or gluons which form 
hadronic jets, and the relevant kinematic variable is the invariant mass rrijj of the two jets thus produced. We assume 
that each tp decays primarily via three-body processes of the form ip — > jjXn- We discuss the alternative possibility in 
which ip decays to the same final state through cascade decays involving a on-shell intermediary in Sect. I VI I Finally, 
we assume here for simplicity that the pair of final-state jets produced via the decay of each particular ip can be 
correctly identified in each event. Such an identification is trivial in situations such as those in which ip is produced 
singly or in which there exist other decay channels for ip with branching fractions similar to ^2 n BR(^ — > jjXn) whose 
decay products are readily distinguishable from jets. However, if no such alternative decay channels exist — i.e., if 
decays of the form ip — > jjXn dominate the width of ip — the analysis becomes significantly more complex because 
of the non-trivial combinatorial issues which arise due to the possibility of incorrect pairings among the final-state 
jets. However, a number of techniques have been developed which can assist in identifying the correct jet pairings in 
such situations. We shall discuss these techniques, and the effect of combinatorial issues in general, in more detail in 
Sect. ED 

We begin our analysis of the invariant-mass distributions which arise in the context of the DDM framework by briefly 
reviewing the characteristics of the corresponding distributions which arise in traditional dark-matter models. This 
will be important for purposes of comparison when we consider the contrasting case of a full DDM ensemble. Towards 
this end, let us consider a traditional dark-matter model in which a parent particle ip likewise decays predominately 
into a three-body final state comprising a dark-matter particle x and a pair of strongly-interacting SM fields (either 
quarks or gluons). The latter appear in the detector as hadronic jets, here labeled j\ and ji. Irrespective of the 
Lorentz or SU(3) C representations of these fields, the differential partial width dT^, = dT{ip — > j\jiX) associated with 
this decay channel is given by the general three-body decay-width formula 

I 2 

^ = 32(2tt)%4 |7W| ^ ' (3 - 1} 

where \M.\ is the matrix element for the parton-level process (averaged over the spin and color states of ip), rrijj is the 
invariant mass of j\ and ji, and mj 2X is the invariant mass of js and x- The normalized invariant- mass distribution 
associated with these decays is therefore given by 

1 dr 4> m ii f "+ 7T-7T 2 ,2 

' \M\dmL y , (3.2) 



dmjj ' 16(27r)3m3 7 m2 _ 1 1 »* 



where the limits of integration are 



m \ + m x _ m % ± \j m % " 2r7i^(m2 + m 2 J3 ) + (ro^ - m 2 ) 



(3-3) 



in the limit that parton masses can be neglected. Note that Eq. (|3.2I) vanishes for rrijj = and for rrijj — — m x . 

The explicit form of the rrijj distribution which arises in any particular model depends on the coupling structure 
between ip, x, and the SM quark or gluon fields. As a concrete example, let us consider the case in which x and ip are 
both fermions, and ip transforms in the octet representation of SU(3) C , while \ is neutral under all SM symmetries. 
In this case, ip can decay via the effective four-fermion interaction term 

£cff = ^fei^Xx^+h.c. (3.4) 

where q denotes a SM quark, t". is the generator of SU (3) in the fundamental representation, the indices i and j 
label the states in the fundamental representation, and a labels the states in the adjoint representation. The spin- 
and color-averaged squared matrix element for an interaction of this form is then given by 

2 

Wf = fjK - ™DK« - O ■ ( 3 - 5 ) 
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FIG. 1: Characteristic invariant-mass distributions in the traditional dark-matter models discussed in the text. The left panel 
shows the normalized dijet invariant-mass distributions associated with the decays of a heavy octet fermion ip with a mass 
my, = 1500 GeV to a quark-antiquark pair and a single dark-matter candidate with a mass m x . The black, red, orange, 
green, blue, and purple curves correspond respectively to m x = {200, 400, 600, 800, 1000, 1200} GeV. The right panel shows the 
distributions for fixed = 1500 GeV and m x = 500 GeV, but for ip and \ with different spins and transformation properties 
under SU(3) C and therefore different coupling structures with the SM fields. The red curve corresponds to the parton-level 
process ip — > qq\ in which ip and x are both fermions, the green curve corresponds to the process ip — > qq\ in which ip and x 
are both scalars, and the blue curve corresponds to the process ip — > gqx in which ip is a fermion and x is a scalar. We observe 
that coupling structure does not induce a dramatic change in the dijet invariant-mass distribution. 



In the left panel of Fig. [TJ we display the rrijj distributions obtained in this case for — 1500 GeV and several 
different choices of m x . Note that each distribution shown features a characteristic mass edge at the kinematic endpoint 
rrijj = m^, — m x , which differs for each m x . However, in most other aspects, the distributions are qualitatively quite 
similar, and in particular have the same overall shape. 

Of course, different coupling structures from the one specified in Eq. (|3.4[) correspond to different functional forms 

for \A4\ , and one might wonder what effect such differences in coupling structure have on the rrijj distribution as a 
whole. In the right panel of Fig. [TJ we compare the rrijj distributions associated with several allowed combinations 
of Lorentz and SU(3) C representations for both parent and daughter particles, which result in different coupling 
structures. The distributions shown correspond to the parton-level process ip — > qq\ for which ip and x are both 
fermions (red), the process ip — > qqx for which ip and x are both scalars (green), and the process ip — > gqx for which 
ip is a fermion and x is a scalar (blue). These decay processes arise in the cases in which ip transforms as an 8, as a 6, 
and as a 15 under SU(3) C , respectively. In each case, we have set m^p = 1500 GeV and m x = 500 GeV. While these 
distributions vary slightly from one another, we once again see that they all have essentially the same shape, which 
includes a characteristic mass edge at rrijj = — m x . Indeed, we observe that the rrijj distributions do not differ 
dramatically among single-particle dark-matter models as a result of differences in coupling structure. 

We now consider how the decay phenomenology of ip in the context of the DDM framework differs from that asso- 
ciated with the traditional dark-matter scenarios discussed above. As we shall see, while each individual constituent 
Xn hi the DDM ensemble contributes to the width of ip and to the overall rrijj distribution in a manner analogous 
to a traditional dark-matter candidate, the collective behavior of these x„ gives rise to distinctly new phenomena. 
For concreteness, we once again focus on the particular case in which ip is a color-octet fermion and the x„ are 
fermions which transform as singlets under the SM gauge group; however, as the right panel of Fig. [TJ attests, the 
shape of the rrijj distributions does not depend sensitively on the spins of the particles involved or the structure of 
the interaction vertex, and the results obtained for different cases will therefore be analogous. We focus on the case 
in which ip — ¥ jjxn decays arise predominately due to a four- fermion interaction analogous to Eq. (|3.4[) : 



-off 



E 



fH^ Q )(x„<Z,) + h, 



(3.6) 



For convenience, we shall henceforth use the symbol 7 to denote the particular scaling exponent 7 Q associated with 
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this operator, so that the operator coefficients c n scale according to the relation 



Cn = C 



too 



(3.7) 



in the manner described in Eq. (|2 . 2|) . We will also assume for simplicity that there are no additional decay channels 
for ip, and that the sum of the branching fractions for all decays of the form ip — ¥ jjxn is effectively unity. Moreover, 
we will assume that the operators which contribute to the decay widths of all Xn with masses m n < to^ are sufficiently 
suppressed so that all such particles amply satisfy the condition in Eq. (|2.4[) and therefore manifest themselves as fx- 
For the coupling structure specified in Eq. (13.6[) . the differential decay width dT^ n = F(ip — > qqXn) for ip decay into 



(3.8) 



any particular Xn is likewise given by Eq. (|3.1[) . and \M\ is given by 



\M\' 



c 2 



)(TO 2 2 „ -TO 2 ) , 



where to J2 „ denotes the invariant mass of \n and j2- Integrating over dm 2 - 2n yields the differential partial decay width 
of ip with respect to rrijj for the particular decay ip — > jjxn- 



dT 



c 2 m 



jj 



yjm*- 2m 2 il {m 2 jj + ml) + (m 



m 



2\ r 



dm 



jj 



96(2tt) 3 to3 A 4 



The partial width for ip —} jjxn is therefore 



F-0n 



384(2tt) 3 A 4 to?, 



to| + to 2 (to 2 j — 2to 2 ) — 2m jj + 



m % ~ 8TO v ,m n (m v , - to„) + 12to^to„ In — f 



m^-TO 2 



jj 



(3.9) 



(3.10) 



Since the products of each such decay mode appear in a collider detector as a pair of jets plus I/Ht, the total differ- 
ential dijet invariant-mass distribution observed is a sum of the mjj distributions for each channel, weighted by the 
corresponding decay branching fraction: 



1 dT 



r^, dm 



jj 



E 

n=0 



i dr 



T^ n dm 



x BR 



(3.11) 



where = T^ n is the total contribution to the decay width of ip from processes of the form ip — ¥ jjXm where 
BR^n = r^n/r^ denotes the branching fraction for the particular decay ip — > jjxn, where n max is the value of n 
corresponding to the heaviest Xn kinematically accessible in ip decay. Note that this result also applies in the case in 
which other decay channels exist for ip with distinguishable final states. The invariant-mass distribution in Eq. p. lip 
therefore depends not only on the parameters Am and 6 (which control the mass spectrum of the Xn and therefore 
the available phase space for each individual decay channel), but also the parameter 7 (which affects the branching 
fractions associated with these channels). 

In order to explicitly illustrate the dependence of the branching fractions BR^ n on 5, 7, and Am, in Fig.[2]we display 
BR^,„ as a function of n. In each of the three panels shown, we have taken m^ = 1500 GeV and too = 200 GeV. The 
curves in the left panel correspond to different choices of Am for fixed 7 = and 8=1; this represents a case in which 
the mass splitting between the Xn is uniform and each of these particles couples to ip with equal strength. In this case, 
the dependence of the branching fraction on n is solely due to the the available phase space for ip — > jjxn decays, 
which decreases monotonically with n up to the kinematic limit — i.e., so long as m n < m^p — above which BR^„ = 0. 
For large values of Am, we see that only a few states are kinematically accessible, and the branching fractions to 
the heavier states are significantly suppressed. By contrast, for small Am, a large number of states are kinematically 
accessible, and BR^,„ decreases quite gradually with increasing n. A similar effect is manifest in the center panel of 
Fig. [H in which the curves in which correspond to different values of 6 for fixed 7 = and Am = 50 GeV. However, 
since increasing S has the effect of increasing the mass gap between the %„, we see that BR^,„ drops more rapidly 
with n when S is large. 

Finally, in the right panel of Fig. [51 we display curves which correspond to different values of 7 for fixed S = 1 and 
Am = 50 GeV. For 7 < 0, we observe that BR^„ decreases monotonically with n as in the left and center panels, 
since such values of 7 simply imply an additional coupling suppression of BR^„ for the heavy states in addition to 
the phase-space suppression discussed above. By contrast, the couplings of the heavy states are enhanced for 7 > 0, 
and the effect of this coupling enhancement and the phase-space suppression compete. As a result, in this case, the 
dominant decay mode for ip may not be to the lightest state in the ensemble, but rather to a more massive state. 
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FIG. 2: The branching fraction BR(i/> — > jjXn), plotted as a function of the DDM-ensemble index n. In all cases we have taken 
my, = 1500 GeV and mo = 200 GeV. In the left panel, we have set 5 = 1 and 7 = 0, and the red, orange, green, blue, and 
purple curves correspond respectively to Am = {600, 500, 400, 150, 20} GeV. In the center panel, we have set Am = 50 GeV and 
7 = 0, and the same colors respectively correspond to 8 — {2, 1.5, 1, 0.75, 0.5}. In the right panel, we have set Am = 50 GeV 
and 5=1, and the same colors respectively correspond to 7 = { — 2, — 1, 0, 1, 2}. Note that the left and center panels show the 
DDM index n on a log scale, while the right panel shows n on a linear scale. 



Having examined how the branching fractions BR^„ depend on the parameters Am, S, and 7 which characterize our 
DDM ensemble, we now turn to examine how the rrijj distributions themselves depend on these parameters. In Fig. [3J 
we illustrate how these distributions depend on 7 and Am for fixed = 1500 GeV, mo = 200 GeV, and 6 = 1. The 
results shown in the left, center, and right panels correspond to 7 = {—1, 0, 1}, respectively. The red, orange, green, 
and blue curves in each panel respectively correspond to the mass splittings Am = {600, 400, 150, 20} GeV, while the 
black curve corresponds to the limiting case in which Am — ¥ 00, or equivalently to the case of a single dark-matter 
candidate with mass m x = mo- 
lt is evident from Fig. [3] that in cases in which 7 < 0, the heavier fields in the DDM ensemble couple more weakly 
to ip than the lighter states, to the result that \ Q dominates in T^. However, as 7 increases, the BRy,„ for n > also 
increase, to the extent that for large 7, decays to the heavier states in the ensemble actually dominate the width of ip. 
In this latter regime, multiple kinematic edges are evident when Am is sizable. While such edges cannot be resolved 
when Am is small, the peak of the distribution nevertheless shifts to smaller values of rrijj. This behavior is due to 
the increased branching fractions of ip to the plethora of heavier \n in the DDM ensemble which are kincmatically 
accessible in ip decays. 

In Fig. [4j we illustrate the dependence of the rrijj distributions on 5 and Am for fixed m^ = 1500 GeV, mo — 
200 GeV, and 7 = 1. The left, center, and right panels respectively correspond to the cases in S — {0.5, 1, 3}. We see 
from the right panel that for large <5, and especially when Am is large, the rrijj distribution is sensitive primarily to 
only the lightest few \n in the DDM ensemble. As a result, the distributions shown are characterized by the presence 
of several identifiable mass edges, each corresponding to one of these light By contrast, we see in the left panel 
that for small 6, a qualitatively different situation emerges: in this regime, the individual mass edges become difficult 
to distinguish, especially for small Am, and the peaks of the invariant-mass distributions shift to lower values of mjj. 

Finally, in Fig. [5l we illustrate the dependence of the mjj distributions on 5 and 7 for fixed = 1500 GeV, 
mo = 200 GeV, and Am = 200 GeV. The left, center, and right panels correspond respective to the cases in which 
5 = {0.5, 1, 2}. As in Fig. 21 for large S we see that the rrijj distribution is sensitive primarily to the lightest states in 
the ensemble and features an identifiable mass edge corresponding to each such state, while for small 6 the mass edges 
become increasingly indistinct. However, we also see from the left panel that the overall shape of the distribution and 
the location of its peak still depend on 7, even in situations in which no individual mass edge can be identified. Thus, 
even in this limit, we find that the shape of the mjj distribution conveys non-trivial information about the structure 
of the DDM ensemble. 

To summarize, we see that DDM models give rise to two characteristic classes of mjj distributions in different 
regimes of model parameter space. Each of these represents a dramatic departure from the mjj distributions typically 
realized in traditional dark-matter models. In the regime in which only a few states contribute significantly to the 
decay width of the parent particle ip, the mjj distribution is characterized by the presence of multiple identifiable 
mass edges. This occurs either when the mass splittings between the heavier states in the ensemble are large, due to 
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FIG. 3: DDM invariant-mass distributions, shown for increasing 7 with fixed = 1500 GeV, mo = 200 GeV, and 5 = 1. The 
results shown in the left, center, and right panels correspond to 7 = { — 1, 0, 1}, respectively. In each panel, the red, orange, 
green, and blue curves correspond to mass splittings Am = {600, 400, 150, 20} GeV, respectively, while the black curve shows 
the result for a traditional dark-matter candidate with m x = mo. 
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FIG. 4: DDM invariant-mass distributions, shown for increasing 8 with fixed = 1500 GeV, mo = 200 GeV, and 7 = 1. The 
results shown in the left, center, and right panels correspond to 8 = {0.5, 1, 3}, respectively. In each panel, the red, orange, 
green, and blue curves correspond to mass splittings Am = {600, 400, 150, 20} GeV, respectively, while the black curve shows 
the result for a traditional dark-matter candidate with m x = mo. 




FIG. 5: DDM invariant-mass distributions, shown for increasing 8 with fixed = 1500 GeV, mo = 200 GeV, and Am = 
200 GeV. From top left to bottom right, the panels shown show results for 8 = {0.5, 1, 2}. In each panel, the red, orange, green, 
blue, and purple curves correspond to 7 = {—2, — 1, 0, 1, 2}, respectively, while the black curve shows the result for a traditional 
dark- matter candidate with m x = mo. 
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either Am or S being sizeable. To some extent, this behavior is not unexpected, but the presence of such identifiable 
mass edges may ultimately provide one clear way of discerning a DDM ensemble experimentally. By contrast, in the 
opposite regime in which a large number of states contribute significantly to the decay width of the parent particle 
-0, the rrijj distribution exhibits no identifiable edges, but instead assumes a unique shape which is markedly different 
from that observed in traditional dark-matter models. Such a shape emerges when mass splittings between the 
heavier states in the ensemble are small, and consequently a larger fraction of that ensemble (consisting primarily of 
the heavier \n) contributes non-trivially to the decay width of ip. 



IV. DISTINGUISHING DDM ENSEMBLES AT THE LHC 



In the previous section, we examined a number of distinctive features which can emerge in the invariant-mass 
distributions associated with ip — > jjxn decays in DDM models. In this section, our primary aim is to assess the 
degree to which the characteristic rrijj distribution associated with a particular DDM model constitutes a distinctive 
signature of non-standard physics in the dark sector — i.e., a signature that cannot be realized in any traditional 
dark-matter model, regardless of the mass m x of the dark-matter candidate. Once again, for concreteness, we focus 
on the operator structure given in Eq. (|3.6|) . In order to determine the statistical significance with which the rrijj 
distribution associated with any particular DDM model (i.e., a particular set of values for the parameters m^,, m , 
Am, 7, 8) is truly distinctive, we compare this distribution to the rrijj distributions associated with a variety of 
different traditional dark-matter models, each with a different value of m x . In particular, we canvass the entire range 

< m x < my, — m jj with a finite step size. Note that we do not similarly scan over all possible spins combinations 
for ip and x> or over coupling structures of these particles to SM quarks and gluons, but rather restrict our analysis 
to traditional dark-matter models in which the parent particle likewise decays via an operator of the form given in 
Eq. (|3.4[) . This is justified because the shape of the rrijj distribution does not depend sensitively on the coupling 
structure, as we have demonstrated in Sect. lIIIl We have nevertheless verified that our results do not differ significantly 
even if alternative coupling structures are incorporated into the analysis. 

The procedure we adopt in comparing any two rrijj distributions is as follows. When rrijj is small, statistics are 
low and residual SM backgrounds are at their largest; we therefore begin by applying a minimum cut on rrijj of the 

form rrijj > m^" , where m^ 1 is some particular invariant-mass threshold. We then partition the range of allowed 

dijet invariant masses m?!j < rrijj < — mo into rn, bins whose widths vary with rrijj. Specifically, the width of 
each bin is taken to be equal to the dijet invariant-mass resolution Arrijj at the minimum rrijj in the bin. We assume 
here that Arrijj is limited predominately by the uncertainty AEj in the measurement of the energies Ej of the jets 
used in reconstructing rrijj, and hence we take (Am jj) /rrijj {AEj)/ Ej. For the ATLAS detector, the jet-energy 
resolution has the rough scaling behavior [ll| 

AEj ( Ej V 1/2 , . 

and the jet-energy resolution at CMS is quite similar [l2| • 

For each value of m x in our survey, we assess the goodness of fit between the rrijj distributions associated with the 
specific DDM model under study and traditional dark-matter models by constructing the % 2 statistic 

xVx) = ^ [Xk - E i m * )] \ (4.2) 

where the index k labels the bin, Xk is the expected population of events in bin k in the DDM model, Ek{m x ) is 
the expected population of events in bin A; in a traditional dark-matter model in which the dark-matter particle has 
mass m x , and a\ is the variance in Xk due to statistical uncertainties. Since the Xk are distributed according to a 
multinomial distribution, it follows that a\ = Xk(l — Xk/N e ), where N e denotes the total number of signal events in 
the sample. Note that the minimum from among these X 2 ( TO x) va lues embodies the degree of discrepancy between 
the rrijj distribution associated with the DDM model and that associated with the traditional dark-matter model 
which provides the best fit to that distribution. We therefore adopt 

Xmin = min{x 2 (m x )} (4.3) 

as our final measure of the distinctiveness of the rrijj distribution associated with the DDM model. We further 
determine the statistical significance associated with a particular value of x min by first comparing this value to a x 2 
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FIG. 6: Contour plots showing the minimum significance level at which the rrijj distribution predicted in the DDM model 
characterized by the parameters my,, mo, Am, 8, and 7 can be differentiated from the rrijj distribution predicted in any 
traditional dark-matter model. In the panels shown, we have fixed my, = 1500 GeV and mo = 200 GeV, and taken m^™ 111 ' = 
200 GeV as our minimum rrijj threshold. The panels in the left column show results for Am = 50 GeV, the center column for 
Am = 200 GeV, and the right column for Am = 500 GeV. In each column, the top panel corresponds to N e = 500, while the 
bottom panel corresponds to N e ~ 1000. 

distribution with n.& — 1 degrees of freedom in order to obtain a p-value. We then take the significance level to be 
the number of standard deviations away from the mean to which the same p-value would correspond for a Gaussian 
distribution. Note that unlike in typical bump-hunting analyses, the p-value in this case is two-sided. 

In Fig. we show how the significance of differentiation varies as a function of S and 7 with Am, too, and 
to^j held fixed. In all of the panels shown, we have set = 1500 GeV and mo = 200 GeV, and we have taken 

TO.j™ ln ^ = 200 GeV as our minimum rrijj threshold and a step size of 100 GeV in our scan over m x . The panels in 
the left column show the results for Ato = 50 GeV, the middle column for Ato = 200 GeV, and the right column 
for Ato = 500 GeV. In each column, the results shown in the top and bottom panels correspond to N e — 500 and 
N e = 1000, respectively. 

The results in the left column indicate that for relatively small Ato, our DDM model can be distinguished from 
traditional dark-matter scenarios at a significance level of 5a or higher for 7 > 0.5 and 6 < 3 with N e > 1000. This is 
due to the fact that the coupling to the heavier Xn are proportionally larger for 7 > 0, and these states increasingly 
dominate the width of "0 as 7 increases. As a result, the peak of the rrijj distribution shifts to smaller smaller values 
of rrijj, and is therefore more readily distinguished from the distribution obtained in traditional dark-matter models. 
By contrast, for 7 < 0, the width of ip is dominated by decays to the lighter x-m and the rrijj distributions obtained 
in this regime resemble much more closely those associated with traditional dark-matter models. This behavior is 
evident in the center and right panels as well. 



12 



We also see from the panels of Fig. [6] that the significance for differentiation depends non-trivially on 6 and Am 
as well. For example, an increase in significance within the 7 < region is obtained for small 5 — an increase 
which is particularly pronounced for small Am. This effect is a consequence of the sheer multiplicity of heavy Xn 
overwhelming the effect of the coupling suppression of ip to these states (which occurs for 7 < 0) and resulting in 
a rrijj distribution which peaks at far smaller values of rrijj than in traditional dark-matter models. Conversely, in 
the opposite regime in which 7 > and 5 is large, a precipitous drop in the significance occurs for 6 > 3 in the 
7 > 0.5 region for Am — 50 GeV, and at even smaller values of 6 for larger Am. This behavior reflects the fact 
that as Am and 5 increase, the mass splittings between the \ n eventually become larger, and consequently fewer and 
fewer Xn remain kinematically accessible in ip decays. In such cases, the contribution from either of these states can 
mimic the m^ distribution obtained for a single dark-matter particle, depending on the branching fractions BR.0O 
and BR^!. Whenever BR^ BR^i or BR^ -C BR^i, only one of these two particles dominates the width of ip, 
and the resulting rrijj distribution resembles that obtained for a traditional dark-matter candidate with m x = mo or 
m x = mi, respectively. However, in situations in which BR^o ~ BR^i and both particles have a salient effect on the 
rrijj distribution, it takes a shape which no single- particle dark-matter candidate can replicate. Since the relationship 
between these branching fractions is primarily governed by 7, there exists a narrow range < 7 < 1 within which the 
significance for differentiation is substantial for large 6 and Am. By contrast, outside of this region, either xo or xi 
overwhelmingly dominates the width of ip, and the significance drops dramatically. Indeed, this behavior is evident 
in the right panels of Fig. [S] 

The plots in each column of Fig. |B] illustrate the effect of N e on the significance of differentiation, and in particular 
on the sensitivity of our results to the total event count. Indeed, we observe that a 5a significance of differentiation is 
achieved over a a substantially larger region of parameter space for N e — 1000 than for N e = 500. Moreover, it turns 
out that the region of parameter space in which such a significance is obtained diminishes rapidly with decreasing 
N e for N e < 500. Thus, we conclude that N e ~ O(500 — 1000) is roughly the size of the data sample for which 
rrijj distributions begin to provide significant resolving power between DDM scenarios and traditional dark-matter 
models. In Sect. [VJ we shall provide an example of a production mechanism which naturally yields such N e values 
for £ int = 30 fb -1 at the y/1 = 14 TeV LHC. 

In Fig. [7J we show how the significance of differentiation varies as a function of Am and 5 with 7, mo, and m^ 

held fixed. The results shown correspond to the same choices of m^, mo, mj™"', and step size for m x as in Fig. [5] 
The left, center, and right panels show the results for 7 = {—1, 0, 1}, respectively. Once again, we see from the left 
panel that for 7 < 0, the density of states as a function of mass within the DDM ensemble must increase quite rapidly 
with n to overcome the corresponding coupling suppression. Both Am and 5 must be extremely small for this to 
occur. As 7 increases, less extreme values for these parameters are required to overcome the coupling suppression, 
as the results shown in the center panel of the figure attest. A horizontal band of slightly increased significance for 
300 GeV < Am < 450 GeV is also evident for 8 > 1.5. This once again corresponds to the region of parameter space 
discussed above, within which xo and Xi are the only states in the ensemble kinematically accessible in ip decays, but 
within which the mjj distribution nevertheless differs significantly from those associated with traditional dark-matter 
candidates because BR^, ~ BR^i. In the right panel, the couplings of the heavier Xn are enhanced relative to those 
of their lighter counterparts, and as a result, the only regions of parameter space shown in which a 5a significance is 
not obtained are those which correspond to the upper right portions of the Am = 50 GeV and Am = 200 GeV panels 
in Fig. |B1 within which only two of the Xn are kinematically accessible in ip decays, and one of these two particles 
overwhelmingly dominates the width of ip. 

Finally, in Fig. [5J we show how the significance of differentiation varies as a function of Am and 7 with 6, mo, and 
held fixed. Once again, the results shown correspond to the same choices of my,, mo, nv™ , and step size for m x 
as in Fig. [51 The left, center, and right panels show the results for 5 = {0.5, 1, 2}, respectively. The results in these 
panels illustrate once again that large 7 generally yields a large significance for differentiation — except in situations 
in which xo and xi are the only Xn sufficiently light to be produced by ip decays, and the condition BR^o ~ BR^,i 
must be satisfied in order for the mjj distribution to differ significantly from that of a traditional dark-matter model 
with either m x — mo or m x — m\. 

Taken together, the results displayed in Figs. [B] through [5] indicate that it is possible to discriminate between a wide 
range of DDM ensembles and more traditional dark-matter candidate on the basis of dijet invariant-mass distributions 
at the LHC. Indeed, we have shown that a statistical significance for discrimination in excess of 5a is obtained for 
N e = 1000 over a broad region of the DDM model parameter space, and regions of parameter space within which 7 
is large or within which Am is small and S < 1 are particularly auspicious. 
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FIG. 7: Contour plots showing the minimum significance level at which the rrijj distribution predicted in the DDM model 
characterized by the parameters m$, mo, Am, 8, and 7 can be differentiated from the mjj distribution predicted in any 
traditional dark-matter model for N e = 1000. The colored regions shown correspond to the same significance intervals as in 
Fig. |6j and the results shown correspond to the same values of m^, mo, and m^ ln \ The left panel shows results for 7 = — 1, 
the center panel for 7 = 0, and the right panel for 7 = 1. 




FIG. 8: Contour plots showing the minimum significance level at which the mjj distribution predicted in the DDM model 
characterized by the parameters m^, mo, Am, 5, and 7 can be differentiated from the mjj distribution predicted in any 
traditional dark-matter model for N e = 1000. The colored regions shown correspond to the same significance intervals as in 
Fig. [6] and the results shown correspond to the same values of m^, mo, and m^ ln \ The left panel shows results for 8 = 0.5, 
the center panel for 5 = 1, and the right panel for S — 2. 



V. DDM PRODUCTION CHANNELS AND EVENT RATES 



Thus far, we have not specified any particular production mechanism for the parent particle ip whose decays give 
rise to the constituent fields of the DDM ensemble. Indeed, all of our results up to this point have been essentially 
independent of the details of the mechanism by which the ip are produced, provided that the SM background can 
be reduced to a negligible level by cuts imposed on the data and that jets associated with the decay of each parent 
particle can be correctly identified. Indeed, they depend only on the total event count N e . This gives our analysis 
a broad range of applicability. We have demonstrated that in a wide variety of situations, iV e ~ O(500 — 1000) is 
sufficient to distinguish a DDM ensemble from a traditional dark-matter candidate at the 5ct significance level. It 
is therefore natural to wonder whether such values of N e might reasonably be expected within the next few years 
of operation at the LHC, and what sorts of parent-particle production processes are capable of yielding such event 
counts. 

In this section, we provide an example of a production process which generically yields event counts of this order 
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FIG. 9: One possible DDM production channel: the pair production of heavy states which decay into four jets and a pair of 
DDM-ensemble constituents. 

for any strongly-interacting parent particle with a TeV-scale mass: the pair-production of ip, as indicated in Fig. |9] 
Indeed, any strongly-interacting exotic can be produced in this manner, via its interactions with the gluon field. 
While the precise value of N e in any particular case depends on the chosen Lorentz and SU(3) C representations of 
ip, on the value of m^, and on the event-selection criteria imposed, we shall now demonstrate that event counts of 
the order assumed in Sect. IIVI can easily be obtained with an integrated luminosity of roug hly £ in t ~ 30 ft) at the 
y/s = U TeV LHC. 

For any process via which ip is produced at colliders, the expected number of signal events N e surviving all cuts 
imposed on the data is given by 



where £j nt is the integrated luminosity, A is the detector acceptance, e S j g is the signal-event-selection efficiency, a^, is 
the cross-section for the production process in question, and is the multiplicity of ip particles produced in each 
event. For pair production, of course, — 2. Note that in assessing e s i g for any particular production mechanism, 
we stipulate that the corresponding event-selection efficiency for the SM background is such that this background is 
negligible. Moreover, in what follows, we shall assume for simplicity that the total branching fraction of ip to final 
states consisting of a pair of jets and one of the Xn is effectively unity. With these assumptions, the contribution to 
N e from pair production depends simply on the total pair-production cross-section a pp ^^ for ip, on the cuts imposed 
on the data (through e S j g ), and on the inherent properties of the detector (through A). 

We begin by evaluating a pp ^^ . The parton- level pair-production cross-sections <jq q ^^ and <j gg ^^ for an exotic 
scalar or fermion which transforms in an arbitrary representation of SU(3) C were calculated to leading order (LO) 
in Ref. [IjJ. Convolving these expressions with the appropriate parton-distribution functions (PDFs) in each case 
yields the LO hadronic cross-section a pp ^^. In this analysis, we use the CTEQ6L1 [14[ PDF set, and impose a 
pseudorapidity cutoff \rj\ < 3 in computing the hadronic cross-sections in order to account for the finite pseudorapidity 
coverage of the LHC detectors. The LO hadronic cross-sections corresponding to several combinations of SU (3) c and 
Lorentz representations for ip for which decays of the form ip — > jjXn can have sizeable branching fractions are shown in 
the left panel of Fig. [TU]as a function of . We see from this figure that in the case in which ip is a Majorana fermion 
transforming in the 8 representation of SU(3) C or a Dirac fermion transforming in the 15 representation, a cross 
section in the range a pp ^^ > 100 fb is attained for all values of < {1100, 1500} GeV at the yfs = 14 TeV LHC. 
While the cross-sections are smaller for the cases in which ip is a scalar field transforming in the 3 or 6 representation 
of SU(3) C , they are still be reasonably large for < 1000 GeV. 

We now turn to address the values of A and e s i g associated with the pair-production of a strongly-interacting 
parent particle. Since e s ; g depends sensitively on the details of the scenario in question and on the particular set 
of event-selection criteria which provides an optimal discriminant between signal and background in that scenario, 
we focus on one representative example case: that in which ip is a color-octet Majorana fermion. For this case, 
representative values for the product ^4e s i g may be taken from results of searches for pair-produced gluinos by the 
ATLAS collaboration [TEI, [l6| in the jets + channel in the regime in which all squarks are taken to be heavy. 
Specifically, we adopt the efficiencies obtained for a number Nj > 3 of high-py jets and rn c ff > 500 GeV, where m e s is 
the scalar sum of fx and the magnitudes pTj of the transverse momenta of the three leading jets in the event, ranked 
by pr- In either case, Ae s [ g lies within the range 0.3 — 0.6 for > 500 GeV. In addition, an estimate of the effects 
of next-to-leading-order (NLO) corrections to a pp ^^ in this case can be obtained using the A-factor formalism, in 
which the LO cross-section is scaled by an overall multiplier. The NLO AT-factors for ip pair production used in this 
analysis were obtained using the Prospino package [13], and vary from K rj 1.7 — 4.4 for 300 < < 2500 GeV at 
y/s = 8 TeV and from K rj 1.5 — 2.1 at yfs = 14 TeV for the same mass range. 

In the right panel of Fig. [TU1 we show the expected number N e of events at the ATLAS detector for the case in 
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FIG. 10: Left panel: The pair-production cross-sections a pp ^^,^ obtained for a variety of relevant combinations of SU(3) C and 
Lorentz representations for a strongly- interacting parent particle ip. The labels shown in the legend indicate the dimension of 
the SU(3) C representation for ip in each case, as well as whether ip is a scalar (as in the 3 and 6 cases shown), a Majorana 
fermion (as in the 8 case), or a Dirac fermion (as in the 15 case). The dashed curves indicate the results for y/s = 8 TeV, while 
the solid curves indicate the results for y/s = 14 TeV. Right panel: the expected number of pair-production events N e (taking 
into account detector acceptance and signal efficiency) for a color-octet Majorana fermion ip, plotted as a function of its mass 
m$ for a variety of different integrated luminosities at the LHC. The red curve corresponds to >Cint = 15 fb _1 at y/s = 8 TeV. 
The two blue curves correspond to & n t = 30 fb _1 (solid curve) and C- m t = 300 fb" 1 (dashed curve) at y/s = 14 TeV. Dotted 
horizontal lines indicating N e = {1, 100, 1000} have also been included for convenience. We thus see that significant values of 
N e can easily be achieved at the LHC through this production mechanism. 

which ip is a color-octet Majorana fermion, plotted as a function of with C- m t = 15 fb _1 for the current LHC 
run at y/s — 8 TeV (red), and with C- mt = {30,300} fb -1 for the projected upcoming run at y/s = 14 TeV (blue). 
These results incorporate the effects of detector acceptance and signal efficiency, as discussed above, and include the 
relevant NLO i^-factors. Dotted lines corresponding to N e = {1,100,1000} have also been included for reference. 
Note that for this choice of representations, is constrained by ATLAS gluino searches [l8[ in the jets+ $t channel 
in the context of supersymmetric models with heavy, decoupled squarks. Such searches likewise place a constraint 
m ip ^ 700 GeV on the mass of any color-octet Majorana fermion which decays to a pair of jets and an invisible 
particle with a branching fraction near unity. 

The results displayed in the right panel of Fig. [TU] demonstrate that the pair-production of strongly-interacting 
parent particles with TeV-scale masses is certainly capable of yielding N e ~ O(500 — 1000) within the first 30 fb _1 
of integrated luminosity at the y/s — 14 TeV LHC — an event count which was shown in Sect. HVI to be sufficient for 
conclusively distinguishing DDM scenarios from more traditional dark-matter models on the basis of mjj distributions. 
Indeed, we see that N e ~ 1000 can be expected for a color-octet fermion parent particle with < 1500 GeV at such 
an integrated luminosity during the planned run at y/s = 14 TeV. Furthermore, even with the integrated luminosity 
Ant ~ 15fb _1 anticipated by the conclusion of the current run at y/s = 8 TeV, an event count of roughly N e ~ 1000 
is obtained for < 800 GeV. While acceptance and signal-efficiency factors will differ somewhat for parent particles 
with different spins and SU(3) C representations, the results shown in the left panel of Fig.[TU]indicate that event counts 
of this order can likewise be expected in a variety of other cases as well. It should be noted, however, that combinatorial 
ambiguities associated with the pairing of final-state jets frequently necessitate the imposition of additional cuts, and 
this can lead to a reduction in signal efficiency. In Sect. I VII we shall discuss various methods through which such 
combinatorial ambiguities can be addressed. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have examined the prospects for observing evidence of dynamical dark matter at colliders, and 
for distinguishing between DDM scenarios and more traditional dark-matter scenarios on the basis of collider data. 
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We have focused primarily on one promising technique: the identification of distinctive features imprinted on the 
kinematic distributions of SM fields produced along with the constituent fields of the DDM ensemble via the decays 
of other, heavier fields which happen to be present in a given DDM model. We have shown that the distributions 
which arise in DDM models can differ substantially from the characteristic distributions which arise in traditional 
dark-matter models in which a single stable, neutral beyond-the-SM particle is responsible for essentially the entirety 
of ficDM- To illustrate the efficacy of our approach, we have examined the prospects for observing a statistically 
significant deviation in the shape of the dijet invariant-mass distributions obtained for a particular representative 
DDM scenario relative to those obtained in traditional dark-matter models. We have demonstrated that throughout 
a substantial region of the parameter space of our scenario, a deviation at the 5cr level is obtained within the first 
30 fb _1 of integrated luminosity at the y/s — 14 TeV LHC. In particular, we have shown that regions of the parameter 
space of this scenario in which 7 > or in which Am is small and 5 < 1 are particularly auspicious for differentiation. 
These findings confirm that the LHC is capable of providing compelling evidence for dynamical dark matter within 
the next few years. 

We emphasize that the quantitative results we have obtained depend primarily on the total production rate for 
the parent particle and not on the particular manner in which this parent particle is produced. Indeed, these results 
apply directly to any situation in which the jets produced by the decay of each parent particle ip in the event 
can be unambiguously identified with that particular ip. In certain cases, such identification is indeed trivial. For 
example, ip may be produced singly or have competing decays. Moreover, even in situations in which tp decays 
predominately via parton-level processes of the form ip — > qqXn, heavy- flavor tagging could assist in the identification 
of the correct pairings in events in which one parent particle decays to a final state involving light quarks and the 
other to a final state involving heavy quarks, e.g., bb. In many situations, however, it is not possible to identify the 
correct pairing of final-state jets based on the properties of the jets alone, and incorrect jet pairings give rise to a 
substantial combinatorial background. Nevertheless, a number of techniques exist for eliminating this combinatorial 
background. These include the hemisphere method [3 US), the use of additional kinematic variables such as mr 2 [2lj |. 
and background-subtraction strategies [22| ■ Indeed, it would be interesting to assess the relative efficacy of these 
techniques in helping to discriminate between DDM models and more traditional dark-matter scenarios in the presence 
of combinatorial ambiguities [23[ . 

We note that we have chosen to focus on the case in which the SM particles produced alongside the Xn by ip decays 
are quarks or gluons solely because large event rates can generically be expected for strongly-interacting ip. However, 
the techniques that we have developed here can also be applied to situations in which ip decays to final states involving 
different combinations of SM particles as well. For example, the prospects for differentiating between DDM scenarios 
and traditional dark-matter models obtained for decay topologies such as ip ~> £ + ^~Xn, where ^ ± denotes a charged 
lepton, will be identical to those obtained for ip — > jjXn, given the same total number of events N e . 

We also note that the techniques outlined here for differentiating between DDM ensembles and more traditional 
dark-matter candidates also have a range of applicability that extends beyond the DDM framework. Indeed, not 
every particle which manifests itself as fj, at colliders necessarily contributes significantly (or at all) to Hcdm, and 
$t signatures of this sort can arise in a number of theories involving multiple neutral particles which are stable on 
collider time scales but not necessarily stable on time scales approaching the age of the universe. As discussed in 
Sect, mi examples of such theories can be found among supersymmetric models with highly compressed spectra [25j . 
scenarios involving universal extra dimensions [8l-fl0l|. and a wide variety of additional extensions of the SM. 

Finally, we note that in this paper we have focused on the case in which ip decays to a final state consisting of 
SM states and one or more of the Xn directly, via an effective multi-body interaction vertex. Alternatively, one 
could also consider the situation in which ip decays to such a final state primarily via decay chains involving on- 
shell intermediate states. Indeed, the kinematic distributions associated with the final-state SM fields in this latter 
case differ significantly from those obtained in the former. For example, in traditional dark-matter models, the rrijj 
distributions associated with ip — > jjx processes which proceed via cascade decays involving an on-shell intermediary 
take a characteristic triangular shape. It would be interesting to examine the corresponding Tfijj distributions which 
arise for DDM ensembles, and whether such dark-matter candidates can likewise conclusively be distinguished from 
traditional dark- matter candidates on the basis of these distributions [2J|. 
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